Dynamic Spin Structure Factor of SrCu 2 (B03)2 at Finite Temperatures 
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Using finite temperature Lanczos technique on finite clusters we calculate dynamical spin structure 
factor of the quasi-two-dimensional dimer spin liquid SrCu2(BOs)2 as a function of wavevector 
and temperature. Unusual temperature dependence of calculated spectra is in agreement with 
inelastic neutron scattering measurements. Normalized peak intensities of the single-triplet peak 
are q-independent, their unusual temperature dependence is analyzed in terms of thermodynamic 
quantities. 

PACS numbers: 75.10.Jm, 75.40.Gb, 75.40.Mg,75.25.+z 



I. INTRODUCTION 

In low-dimensional quantum spin systems quantum 
fluctuations often lead to disordered ground states 
that exhibit no magnetic ordering and a gapped, non- 
degenerate singlet ground state. Such states, also called 
spin liquids, are realized in one dimension in dimerised 
or frustrated spin chains, even-leg spin ladders as well 
as in the two dimensional Shastry-Sutherland (SHS) 
model, 1 SrCu2(B03)2 is a quasi-two-dimensional spin 
system with a unique spin-rotation invariant exchange 
topology that leads to a singlet dimer ground stated 
Since this compound represents the only known realiza- 
tion of the SHS model, it recently became a focal point of 
theoretical as well as experimental investigations in the 
field of frustrated spin systems. Consequently, many fas- 
cinating physical properties of SrCu2(B03)2 have been 
discovered. Increasing external magnetic field leads to 
a formation of magnetization plateaus^ which are a 
consequence of repulsive interaction between almost lo- 
calized triplets. Weak anisotropic spin interactions can 
have disproportionately strong effect in a frustrated sys- 
tem. It has recently been shown, that the inclusion of 
the nearest neighbor (NN) and next-nearest neighbor 
(NNN) Dzyaloshinsky-Moriya (DM) interactions is re- 
quired to explain some qualitative features of the spe- 
cific heat near the transition from the spin dimer to the 
spin-triplet state, as well as to explain the low frequency 
lines observed in electron spin resonance experiments in 
SrCu2(BO 3 )2. 5 i 6 ' 7 ' 8 ' 9 ' 10 i 11 

The existence of the spin gap, almost localized spin- 
triplet excited states, as well as the proximity of a 
spin- liquid ground state of SrCu2(BC>3)2 to the or- 
dered antiferromagnetic state, lead to rather unusual 
low-temperature properties emerging in inelastic neutron 
scattering (INS)ji2*i£ Raman scattering (RS)^ as well as 
in electron spin resonance (ESR) experiments. 7 In partic- 
ular, INS normalized peak intensities of single-, double- 
and possibly triple-modes show a rapid decrease with 
temperature around 13 K, well below the value of the spin 
gap energy A ~ 34 K. In addition, authors of Ref>i£ show, 
that properly normalized complement of static uniform 



spin susceptibility, obtained with almost identical model 
parameters as in the present work^ nearly perfectly fits 
their experimental data. Similar behavior is found in RS 
data where a dramatic decrease of Raman modes, rep- 
resenting transitions between the ground state and ex- 
cited singlets, with increasing temperature at T <C A is 
observed. Moreover, at T ~ A all RS modes become 
strongly overdamped^ 

Numerical simulations of dynamical spin structure fac- 
tor based on exact diagonalization on small clusters at 
zero temperature show good agreement with INS data^ 
Recently developed zero-temperature method based 
on perturbative continuous unitary transformationsi&i 7 - 
gives very reliable results for the dynamical spin struc- 
ture factor of the SHS model since the method does not 
suffer from finite-size effects. The method is, however, 
limited to calculations at zero temperature and, at least 
at the present stage, it does not allow for the inclusion 
of DM terms. 

The aim of this work is to investigate finite temper- 
ature properties of the dynamical spin structure factor 
of the SHS model using the finite temperature Lanczos 
method (FTLM)^SiiS and to compare results with INS 
dataJ2il2 In our search for deeper physical understanding 
of spectral properties of the SHS model at finite tempera- 
tures we compare those with thermodynamic quantities, 
such as: the specific heat, entropy, and uniform static 
magnetic susceptibility, which we further compare with 
analytical results of the isolated dimer (DIM) model. We 
finally present results of the q-dependent static magnetic 
susceptibility as a function of T . 



To describe 
SrCu 2 (B0 3 )2 



II. MODEL 

the low-temperature properties of 
we consider the following Heisenberg 
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Hamiltonian defined on a 2D Shastry-Sutherland lattices^ 
H s = J^Si-Sj + J'^Si-Sj 

<u> <u>' 

+ £ D'-(SiXSj). (1) 

(i-j)' 

Here, and (i, j)' indicate that i and j are NN and 
NNN, respectively. A recently published high-resolution 
INS measurments on SrCu 2 (B03)2i^ motivated us to 
choose J = 76.8 K (in units of ks) and J'/J = 0.62. 
A slightly different choice of parameters (J = 74.0 K, 
J'/J = 0.62) has been used previously in describing 
specific heat measurements 8 - and ESR experiments on 
SrCu2(B03)2>^i We should note, however, that this fine- 
tunning of parameters leads to effects, visible only on 
small energy scales, thus leaving previous calculations^^ 
practically unaffected. In addition to the Shastry- 
Sutherland Hamiltonian, H s includes DM interactions to 
NNN with the corresponding DM vector D'. The arrows 
indicate that bonds have a particular orientation as de- 
scribed in Ref£. Its value, D' = 1.77 K z, successfully ex- 
plains the splitting between the two single-triplet excita- 
tions observed with ESR^i and INS measurements ,A2iiii 

As pointed out in Refs&ii, a finite NN DM term 
should also be taken into account to explain specific heat 
data and ESR experiments. We have chosen to omit this 
term since it does not significantly affect results of the 
dynamical spin structure factor at a non-zero value of 
the wavevector. We have chosen the quantization axis z 
to be parallel to the c-axis and x to the a-axis pointing 
along the centers of neighboring parallel dimers. 

We use the FTLM based on the Lanczos procedure of 
exact diagonalization, combined with random sampling 
over initial wave functions. For a detailed explaination 
of the method and a definition of method parameters 
see Refsi 18 ' 19 . All the results are computed on a tilted 
square lattice of N — 20 sites with Mi = 100 first and 
M2 = 250 second Lanczos steps, respectively. The full 
trace summation over N st — 2 N states is replaced by 
a much smaller set of R ~ 500 random states giving the 
sampling ratio R/N st ~ 5 • 10~ 4 . Comparing FTLM with 
the conventional Quantum Monte Carlo (QMC) methods 
we emphasize the following advantages: (a) the absence 
of the minus-sign problem that usually hinders QMC cal- 
culations of frustrated spin systems, (b) the method con- 
nects the high- and low-temperature regimes in a contin- 
uous fashion, and (c) dynamic properties can be calcu- 
lated straightforwardly in the real time in contrast to em- 
ploying the analytical continuation from the imaginary 
time, necessary when using QMC calculations. Among 
the shortcomings of FTLM is its limitation to small lat- 
tices that leads to the appearance of finite-size effects 
as the temperature is lowered below a certain T < Tf s . 
Due to the existence of the gap in the excitation spec- 
trum and the almost localized nature of the lowest triplet 
excitation, we estimate Tf s ~ 1 K at least for calcu- 
lation of thermodynamic properties^ Finite-size effects 



also affect dynamical properties as, e.g., the dynamical 
spin structure factor, which is (even at finite tempera- 
tures) represented as a set of delta functions. In partic- 
ular, finite-size effects affect the frequency resolution at 
low temperatures, while at higher temperatures as more 
states contribute to the spectra, its shape becomes less 
size dependent. 

We should stress that FTLM was in the past success- 
fully used in obtaining thermodynamic as well as dy- 
namic properties of different correlated models as are: 
the t-J model; 18 ' 19 the Hubbard model^S as well as the 
SHS modeU 



III. NUMERICAL RESULTS AND 
COMPARISON WITH EXPERIMENT 

A. Dynamical Spin Structure Factor 

For comparison with the INS data, we compute the 
dynamical spin structure factor for 11 — x,y, z 

00 

Snfau) = ReJdt e- 4 (S$(t)^ q (Q)) ; 


5 q = ^E^( R » + r «) e4q(R!+r ° ) ' ( 2 ) 

where i runs over all unit cells of the lattice and r a spans 
four vectors forming the basis of the unit cell that con- 
tains two orthogonal dimers . For the details describing 
interatomic distances we refer the reader to Refii&. The 
average in Eq. © represents the thermodynamic average 
which is computed using FTLMii 8 . 

In Fig. ^a) we first present the spin structure fac- 
tor S zz (q,ui) for different temperatures. We should note 
that due to a finite value D' z = 1.77 K spin rotational 
invariancc of the Hamiltonian, Eq. is broken, i.e. 
S xx (c[,uj) = S yy (q,uj) ^ S zz (q,uj). Since D' z < J, J', 
the effect of broken symmetry is, at least within our nu- 
merical precision, negligible for energy resolutions much 
larger than the value of anisotropic interaction, Ato 3> 
D' z , yielding nearly identical results for the three compo- 
nents ll — x,y, and z of ^^(q, u>). Since the spectra con- 
sist of a set of delta functions, we have artificially broad- 
ened the peaks with a Gaussian form with a = 0.4 meV, 
to achieve the best fit with INS measurements^ Two 
peaks (I and II) are clearly visible at low temperatures 
T=2K<A around to ~ 3 meV and 5 meV, associated 
with transitions to single- and double-triplet states. A 
broader peak (III) around lu = 9 meV represents tran- 
sitions to triple-triplet states. Results at T = 2 K are 
consistent with previous T = simulations4£*i& Increas- 
ing temperature has a pronounced effect on S zz (q, ui), 
manifesting in a rapid decrease with temperature, at tem- 
peratures even far below the value of the gap. Quantita- 
tively, at T = 24 K the peak structure almost completely 
disappears. In Fig. ^b) we present comparison of our 
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Figure 1: (Color online) (a) Spin structure factor S Z z(<l,ui) for 
q = (2, 0) where q is given in units of the reciprocal lattice 
vectors— vs. to for different values of T, and (b) comparison 
with INS measurements from Refi-. Transitions to single-, 
double-, and triple-triplet states are indicated with Roman 
numerals, (c) Sum of the transverse S yy (q,Lj) and longitudi- 
nal 5 zz (q, ui) contributions of a single-triplet mode compared 
to the high-resolution data for q = (—2,0) and T = 1.4 K 
from RefA We use units where h = 1. 



numerical data scaled and shifted along the vertical axis 
to compensate for experimental background for two dif- 
ferent temperatures along with experimental values from 
Refit 

In order to investigate magnetic anisotropy effects one 
has to turn to high-energy resolution calculations with 
frequency precision comparable to the magnitude of the 
DM interaction. On this frequency scale we expect to find 
substancial difference between longitudinal and trans- 
verse components of the spin structure factor. For this 
reason we have included the transverse component of the 
spin structure factor according to the relation for the dif- 
ferential cross section 



d 2 a 



dfl dui 

\1U 



OC 



(3) 



For a given direction of the neutron momentum trans- 
fer, e.g., q = (—2,0) as used to obtain high-resolution 
INS data presented in Fig.^c), Eq. © reduces to a sum 



of the transverse and longitudinal part, d 2 a / dSl duj oc 
S yy (c[,uj) + S zz (c[,uj). In Fig. QJc) both contributions 
as well as their sum (Eq. (PJ)) are plotted against the 
high-resolution INS datai—i The best fit is achieved for 
J = 76.8 K, J' = 47.6 K, D' z = 1.77 K, and artificial 
broadening of the Gaussian form with a = 0.05 meV. The 
splitting A ~ 0.32 meV between the outer two modes be- 
longing to S 2 = ±1 single-triplet excitations originates 
from a finite value D' z yielding 4D' Z /A ~ 1.9 for the 
renormalization of the bandwidth due to quantum fluc- 
tuations. This is as well in agreement with the estimate 
of Cepas et al.A 

For comparison we also present S zz (q, ui) of the sim- 
plistic DIM model with J = A = 34 K, where J is chosen 
in such a way that SHS and DIM model share identical 
energy gaps between the singlet ground state and the ex- 
cited triplet state. In the latter case analytical expression 
for S zz (q, ui) can be straightforwardly derived 



S zz (q,uj) = vrA(q) (S(u - J) + e" 13 ' 1 5{u 
+ 2nB((i)e- pj S(uj) 

with A(q) and B(q) given by 



J)) 



(4) 



a - sin2 ^ -%)+ sin2 rtv* + 1v) (K\ 



B(q) = 



4(1 + Ze-P J ) 
cos 2 i](q x - q y ) + cos 2 -q (q x + q y ) 
4(1 + 3e-/") 



(0) 



and rj = 0.717. [Note also that in the limit T — * oo 
JZo S zz (q., u>) duj = tt/4.] At T = S zz (q., u>) consists of 
a single delta function at ui — J weighted by A(q). 12 This 
peak corresponds to peak I in the SHS model, while peaks 
II and III do not have their counterparts in S' 2Z (q, ui) of 
the simplistic DIM model. With increasing T peaks at 
u) = — J and u) = appear, weighted by 7rA(q) exp(— j3J) 
and 27ri?(q) exp(— (3 J), respectively. 

In Fig. |21 we present a map of S zz (q, ui) for different 
temperatures. Peak I with almost no dispersion is clearly 
visible at u> ~ 3 meV with its highest intensity located 
near q x ~ 2. Peak II, located at u> ~ 5 meV is also visi- 
ble and similarly shows little dispersion. Its intensity is 
as well maximal near q x ~ 2. Note that due to a small 
system size 5' zz (q, ui) plots were calculated at only a few 
discrete values of q x , i.e., q x = 0.0, 1.0,2.0 and 3.0. The 
geometry of the tilted square lattice with N = 20 sites 
excludes half- integer values of q x . This fact prevents us 
to directly compare our intensity plot results for the spin 
structure factor with the ones shown in Refi^, where the 
dispersion of the lowest triplet mode, attributed mainly 
to the transverse part of the spin structure factor, is 
clearly seen. We have calculated the transverse com- 
ponent Sy y (q,u>) but since it does not differ consider- 
ably from S zz (q, ui) on a given energy scale we do not 
present it in Fig. A final map, shown in Fig. [2 was 
obtained by interpolation between given integer values of 
q x points. These results are roughly consistent with mea- 
surements by Gaulin et al.A With increasing T peaks I 
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Figure 2: (Color online) Intensity plot of the spin structure 
factor Szz (q, w) in the q x — u> plane for different values of T 
as indicated in the figures. The intensity scale is logarithmic 
and a — 0.2 meV is used to broaden the spectra as a function 
of frequency ui. Note that due to the small system size only 
integer values of q x were numerically accesible. 



and II rapidly decrease (more quantitative analysis of the 
temperature dependence follows in the next subsection), 
while visible response due to elastic transitions among 
identical multiplets starts developing around to — and 
q x = 0. 



B. Normalized Peak Intensities 

With the purpose to further quantify the agreement of 
our calculations with the experiment we present in Fig. [3] 
the normalized peak intensities I/Io of the two peaks (I 
and II) as functions of temperature along with the mea- 
sured data taken from Refii 3 -. To avoid contributions 
from the background at higher temperatures, peak in- 
tensities were measured from their values at T — 20 K, 



/0 S 2Z (q,^ P ,0)-S zz (q,u;p,20K)' 



(7) 



with uip = 3.0 meV and uip = 5.0 meV for peak I 
and II, respectively. Gaussian broadening with a = 0.4 
meV was used to obtain peak values of 5 Z2 (q, uip, T). 
A similar temperature behavior is observed as in INS 
measurements manifesting itself in a rapid decrease 
of both peak intensities with temperature for T far be- 
low the gap value A = 34 K. Taking this fact into ac- 



count, the agreement between experimental values and 
numerical calculations of I /Iq is reasonable even though 
not ideal. However, as already suggested by Gaulin et 
al.^ nearly perfect agreement between experiment and 
rescaled complement of the uniform static susceptibil- 
ity 1 - x = 1 - Xo(T)/xo(T = 20 K) is found where 
Xo = (St ot 2 ) /NT, and S% ot represents the z-component 
of the total spin. We were unable to find a direct ana- 
lytical connection between the two quantities, i.e., 1 — \ 
and I/Io, on the other hand the difference between the 
above mentioned quantities in numerical results are obvi- 
ous from Fig. [3J Furthermore, analytical calculations of 
1 — x an d I/Io ° n DIM model also indicate a different T 
dependence. This leads us to the conclusion, that nearly 
perfect agreement with 1 — x and experimental results of 
Ref. 13 may be accidental. 

In the inset of Fig. [3] we show the integrated in- 
tensities of the two peaks defined as I(oj\,u>2) / Io — 
duiS zz (q, u))/{SqS-q) where the limits of integration 
denning integrated peak intensities are defined as follows: 
h = 1(2 meV,4 meV) and I u = 7(4 meV, 6 meV) for 
peaks I and II, respectively. We observe a distinctive 
difference in temperature behavior between // / Iq on the 
one hand and In/Io on the other. While Ii/Iq substan- 
tially decreases with increasing temperature similarly to 
I /Iq, indicating on a considerable shift of the spectral 
weight away from transition I, Ih/Iq shows even a slight 
temperature increase. We suggest that this difference is 
caused by a different nature of the transition from the 
ground state to the localized triplet (peak I) in contrast 
to transitions to states near or else within continuum 
(peak II) . This behavior is as well in agreement with INS 
measurements— that show peak I being only resolution 
limited while peaks II and III show intrinsic lincwidths. 
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Figure 3: (Color online) Normalized peak intensities I/Io of 
peaks I and II vs. T from Refpi^ (filled symbols), FTLM 
results (lines). The inset: relative integrated intensities I/Io 
of peaks I and II. 

We now explore the q-dependence of peak intensities. 
In Fig. 01 we present normalized values of peak intensities 
vs. T for various values of q x at fixed q y — and find 
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nearly perfect scaling of I/Iq for peak I, Fig. Ufa), while 
scaling breaks down at q x = 3.0 for peak II, Fig. 0Jb). 
Such behavior is characteristic also for the simpler DIM 
model that possesses a single temperature scale J. This 
result suggests that a single temperature scale is respon- 
sible for the T dependence of peak I for all different values 
of q. In the insets of Figs. Ola) and^Jb) we present ab- 
solute values of peak intensities for different values of q x . 
Intensities of both peaks I and II reach their maximum 
values at low-T near q x ~ 2.0. Taking into account our 
rather poor resolution in the g^-space, we find these re- 
sults to be roughly consistent with recent high-resolution 
INS measurements by Gaulin et 
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Figure 4: (Color online) Normalized peak intensities I /la of 
peaks (a) I and (b) II vs. T calculated at different values 
of q = (q x , 0). In insets unrenormalized peak intensities mea- 
sured from the peak intensity position at T = 20 K are shown. 



also observe a rapid drop of I/Iq [see Fig. [3], furthermore, 
this temperature also coincides with the steepest ascent 
of s. The peak in c, located well below the value of the 
gap, T max <~ 0.24A, is a consequence of excitations from 
the ground state to localized singlet and triplet states. 
This peak is followed by a broad shoulder above T > 15 K 
which is due to excitations in the continuum. Note that 
c, obtained using the same method and for slightly dif- 
ferent values of J, J' and D, apart for additional DM 
terms, fits measured specific heat data of SrCu2(BOa)2 
in a wide range of applied external magnetic fields^ For 
comparison we also present c and s of the DIM model 
with J = A = 34 K that can be solved analytically. We 
should point out that even in a simple DIM model c peaks 
well below the gap value, i.e., T max = 0.35A ~ 11.9 K. 

In Fig. Efb) we present the uniform static spin sus- 
ceptibility, xo = (Stot 2 )/NT 7 where S* ot represents the 
z-component of the total spin. While comparison of xo 
with experimental data was presented elsewhere^ in this 
work we present it along with other thermodynamic prop- 
erties just to gain a more complete physical picture of the 
temperature dependence of 5 zz (q, to). The steepest in- 
crease in xo vs. T coincides with the peak in c and, at 
least approximately, with the steepest decreases of I/Iq, 
presented in Fig. |3J Due to the existence of the spin 
gap, Xo(T < 5 K) ~ 0, where the temperature interval 
T < 5 K in turn corresponds to the plateau of I/I Q seen 
in experimental results of peaks I and II as well as in our 
numerical simulations. 

We would like to make some general remarks on com- 
paring thermodynamic properties of the SHS model and 
the DIM model with identical gaps between the ground 
state and first excited states. Such a direct comparison 
may assist in understanding the influence of spin frustra- 
tion and the proximity of gapless excitations in the SHS 
model on its thermodynamic properties. In particular, 
the specific heat c of the SHS model peaks at lower tem- 
perature than c of the DIM model and shows two max- 
ima in contrast to a single, Schottky-like maximum seen 
in DIM model. The entropy of SHS model reveals slower 
increase with T than that of the DIM model. And finally, 
the peak value of xo is almost three times lower in the 
SHS model than in the DIM model which in turn implies 
that spin fluctuations, (S z2 ), [see the inset of Fig.[SJb)] 
of the SHS model are suppressed in comparison to the 
DIM model. 



C. Thermodynamic Properties 

Next, we will connect the spectral data with thermo- 
dynamic properties of the Hamiltonian defined in Eq. Q. 
For this reason we present in Fig. OJa) the specific heat 
(per spin) c = T(ds/dT) = ((H 2 ) - (H) 2 )/NT 2 and the 
entropy density s — lnZ/N+(H) /NT, where Z is the sta- 
tistical sum and N is the number of spins in the system. 
Specific heat peaks around T = T max ~ 8 K, where we 



D. Static Spin Susceptibility x(q) 

Finally, we present in Fig. [5] the static spin suscepti- 
bility 



oc 

x(q) = \v J 



x" (q, 



du>, 



x"(q,fcO = (l-e-**)S„(cbw), 



(8) 



(9) 
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Figure 5: (Color online) Thermodynamic properties of the 
SHS model (full lines) and DIM model with J = 34 K (dashed 
lines): (a) specific heat c (entropy s in the inset) vs. T and 
(b) uniform static susceptibility xo vs. T. For completeness 
and to facilitate comparison we present in (c) normalized peak 
intensity of peak I I/Io vs. T measured from values at T = 
20 K as shown in Fig. |3] 

as a function of q = (q x , 0). Besides fulfilling theoretical 
interest, x(q) can also be used to compute, e.g., spin- 
spin nuclear relaxation rate I/T2. Along x(q) of the 
SHS model we present for comparison results for the DIM 
model where analytical result can be readily obtained us- 
ing Eqs. (QJ, ©, and ©: 

X (q) - 2A(q) (1 - e-^) /J + 2B( q )e^ J p (10) 

where A(q) and B(q) are defined in Eqs. (0 and 0. At 
low-T, i.e., T < 5 K, x{°L) vs - ^ is nearly T-independent 
which is a consequence of the spin gap. As a func- 
tion of q x it reaches its maximum value near q x ~ 2.0 
in accord with the prediction of the DIM model result, 
Eq. (|10fl . Observed T-dependence is again similar to 
the DIM model prediction. At higher temperatures, i.e., 
T > 300 K, x(l) merges with universal, q-independent 
form, i.e., x(q) — 1/4T. 

IV. CONCLUSIONS 

In conclusion, we have computed dynamical spin struc- 
ture factor at finite temperatures. Frequency dependence 
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Figure 6: (Color online) x(q = (<lx, 0)) vs. T for several values 
of q x . Circles connected with lines present results of the SHS 
model, lines present results of the DIM model, Eq. (|I0I with 
J = 34 K. 

of S zz (q, uj) at T = 2 and 24 K agree reasonably well 
with INS measurements^ on a large energy scale de- 
spite rather poor frequency resolution caused by finite- 
size effects. High-resolution data for the lowest triplet 
excitation—, on the other hand, is almost perfectly cap- 
tured by calculated transverse and longitudinal compo- 
nents, S yy (q, uj) + S zz (q,u>), showing the influence of 
anisotropy present in the system. Comparison of results 
obtained on systems with = 16 and 20 sites reveals 
that the positions of peaks I and II are reasonably well 
reproduced on N — 20 system, while the peak III posi- 
tion and its structure are less accurate. We should also 
note that although in a 16-site system finite-size effects 
are somewhat more pronounced, this cluster enables one 
to calculate half-integer values of q x as well. By fitting 
the data from Refii^ for all measured q x on a system 
with N = 16 sites, we have concluded that additional 
NN in-plane as well as 'forbiden' NN out-of-plane DM 
interactions are required to succesfully explain the struc- 
ture and dispersion of peak I. However, the best fit is 
obtained with slightly different set of parameters as used 
in this work. 

Temperature dependence of normalized peak intensi- 
ties I/Io agrees well with INS measurements Our 
calculations predict that I/Io of peak I should be q- 
independent. Such behavior is in agreement with the 
DIM model prediction for peak I, while peak II is any- 
how absent in this simplistic model. Our results are 
thus consistent with a proposition that a single temper- 
ature scale is responsible for the T dependence of peak 
I for all different values of q. This statement does not 
take into account a possible small dispersion of peak I 
due to DM interaction or (and) due to high-order pro- 
cesses in J' /J. 1 - From comparison of temperature de- 
pendence of I/Iq with thermodynamic properties it is 
obvious that strong T-dependence of I/Io, occurring well 
below the value of the spin gap, is in accord with strong 
T-dependence of other thermodynamic properties. The 
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temperature of the steepest decrease of I/Iq coincides 
with the peak in c and the steepest increase of s as well 
as of xo- 

There is obviously a need for further, less finite-size de- 
pendent calculations that will clarify many unanswered 
questions as are, e.g., the role of DM terms in explain- 
ing small dispersion of peaks I and II observed in high- 
resolution INS experiments^, an explanation of unusual 
temperature dependence of ESR linesi that seem to de- 
crease in width as the temperature increases, the occur- 
rence of magnetization plateaus, etc. Nevertheless, the 
main features of temperature-dependent dynamic prop- 



erties of the SHS model seem to be well captured by the 
FTLM on small lattices which is in turn reflected in a 
good agreement with experiments. 
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